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1  Introduction 


Background 

Construction  materials  account  for  a  major  portion  of  building  costs,  affecting  the  costs 
of  initial  construction,  continued  maintenance,  habitability,  and  eventual  demolition 
of  a  building.  The  initial  selection  of  a  construction  material  may  depend  on  a  number 
of  complex  and  often  intangible  factors,  but  the  total  initial  and  long-term  costs  of 
using  any  building  material  system  is  one  of  the  most  important  parameters.  Materi¬ 
als  with  the  lowest  initial  cost  often  result  in  such  high  maintenance  costs  that  any 
initial  cost  savings  can  be  erased  much  earlier  than  expected.  A  key  to  cost  reduction 
through  materials  selection  is  knowledge  of  the  comparative  durability  of  various 
building  materials,  and  how  their  durability  affects  the  total  building’s  total  life-cycle 
cost. 

The  service  life  of  a  structure  is  usually  considered  during  design  only  through  the 
specification  of  the  return  period  of  the  most  extreme  environment  the  structure  is 
expected  to  endure.  For  example,  a  structure  may  be  designed  to  withstand  a  100-year 
earthquake  or  wind  storm.  The  selection  of  materials  t3q)ically  is  based  on  the  notion 
that,  with  some  safety  factor  built  in,  the  structure  will  resist  the  extreme  event  in  a 
virgin  state.  A  good  engineer  considers  the  implications  of  material  durability  and 
maintenance  in  the  design  of  a  structure  even  though  current  Army  design  require¬ 
ments  do  not  explicitly  insist  on  it.  The  current  state  of  the  art  in  engineering  design 
for  durability  is  based  on  rather  coarse  decision  points  such  as  “steel  corrodes  while 
plastic  does  not,”  or  “concrete  scales  while  steel  does  not.”  Such  anecdotal  knowledge 
does  not  really  help  the  engineer  evaluate  the  tradeoffs  in  materials  selections  because 
it  is  too  superficial.  The  significance  of  materials  selection  on  overall  Army  operational 
costs  can  be  seen  in  corrosion  cost  figures,  for  example,  which  exceed  $300  million 
annually  (Hahin  1977).  Corrosion  costs  control  are  often  most  easily  reduced  during 
facility  design  through  the  selection  of  alternative  materials. 

The  lack  of  long-term  experience  with  promising  new  materials  tends  to  inhibit 
application  of  new  materials  that  might  significantly  reduce  long-term  facility  costs. 
Facility  designers  and  specifiers  need  a  model  that  can  predict  the  durability  of 
various  building  materials  as  used  for  a  specific  application  at  a  specific  location.  An 
effective  model  of  this  t5q)e  also  would  analyze  what  role  or  priority  materials 
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durability  should  play  in  the  selection  process.  Such  a  model,  when  used  during 
concept  design  in  a  concurrent  engineering  collaborative  decision  environment,  would 
enable  engineers  to  accurately  and  efficiently  compare  the  life-cycle  costs  of  conven¬ 
tional  materials  to  alternative  ones.  The  model  could  help  the  Army  maintain  the  life 
expectancy  of  structures  at  a  lower  life-cycle  cost. 

The  U.S.  Army  Construction  Engineering  Research  Laboratories  was  tasked  by 
Headquarters,  U.S.  Army  Corps  of  Engineers,  to  develop  a  Building  Materials 
Durability  Model  (BMDM). 


Objective 

The  objective  of  this  work  was  to  investigate  concepts  for  optimal  selection  of  building 
materials,  for  possible  incorporation  into  a  decision-support  tool  for  facility  designers 
using  an  agent-based  collaborative  design  network. 


Approach 

The  researchers  investigated  modeling  techniques  that  consider  predictive,  long-term 
materials  performance  in  service,  maintenance  and  repair  (M&R)  life-cycle  implica¬ 
tions,  and  differing  construction  environments. 

A  probabalistic  modeling  approach  was  used.  Damage  models  were  developed  for 
several  common  construction  material  systems,  both  for  structiiral  elements  and 
cladding  elements.  The  model’s  algorithms  were  embedded  into  a  spreadsheet  soft¬ 
ware  program  for  desktop  computer.  Token  building  elements  were  subjected  to 
differing  environments  for  a  number  of  iterations  in  a  computer-driven  Monte  Carlo 
simulation,  and  material  performance  was  evaluated.  Repair  and  failure  criteria  were 
evaluated  after  each  iteration,  and  the  accumulated  cost  of  durabiUty  calculated.  A 
presentation  format  for  durability  comparison  of  different  construction  materials  was 
developed.  Finally,  BMDM  was  applied  to  four  materials — two  structural  materials 
and  two  cladding  materials.  Durability  costs,  in  terms  of  net  present  value,  were 
caluclated  for  two  different  natural  environments. 


Mode  of  Technology  Transfer 

This  work  will  feed  into  the  Architect’s  Associate  program,  USACERL’s  concurrent 
engineering  research  initiative.  The  Architect’s  Associate  technology  transfer  plan  will 
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include  the  formation  of  cooperative  research  and  development  agreements  (CRaDAs) 
for  transfer  of  applications  and  development  of  an  agent  collaboration  language.  The 
agent  collaboration  language  will  be  proposed  for  inclusion  by  ANSI  and  ISO*  in  the 
Product  Data  Exchange  Specification  (PDES)  and  the  Standard  for  the  Exchange  of 
Product  Model  Data  (STEP). 


Metric  Conversion  Factors 

U.S.  standard  units  of  measure  are  used  throughout  this  report.  A  table  of  metric 
conversion  factors  is  presented  below. 


1  in. 

=  25.4  mm 

1  sq  in. 

=  645.16  mm® 

1  lb 

=  3.78  L 

°F 

=  (°Cx1.8)+32 

*  ANSI  is  the  American  National  Standards  Institute;  ISO  is  the  International  Standards  Organization. 
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2  Description  and  Deveiopment  of  BMDM 


Durability  Model  Concept 

Suppose  that  it  is  possible  to  simulate  the  evolution  of  all  environmental  stimuli  that 
are  important  to  the  durability  of  a  building  component.  Suppose  further  that  it  is 
possible  to  analytically  model  the  damage  processes  that  can  occur  in  the  building 
component  and  hypothesize  various  repair  strategies.  Given  these  assumptions,  it 
would  be  possible  to  simulate  the  lifetime  of  that  building  component  over  time  within 
the  framework  described  in  Figure  1*.  In  this  framework  the  environmental  stimuli 
T(t),  d(t),  and  c(t)  are  specified  as  explicit  (possibly  random)  functions  of  time.  The 
component  includes  a  number  of  descriptors,  including  its  material  properties  and  its 
geometric  configuration.  The  quantities  x  represent  a  collection  of  damage  variables 
recording  loss  of  material,  change  in  material  property,  loss  of  aesthetic  appeal,  and 
other  factors  that  evolve  according  to  the  damage  rate  functions  f.  With  an  idea  of  the 
component’s  initieil  cost,  the  criterion  y  can  be  established  to  indicate  failure  or 
to  trigger  repair  and  establish  the  means  to  specify  a  repair  strategy  and  accumulate 
durability  cost  (i.e.,  initial  cost  plus  the  sum  of  the  repair  costs). 

During  each  time  period,  the  amount  by  which  the  material  service  properties  have 
degraded  is  determined  according  to  the  damage  model.  The  degraded  building 
element  is  checked  for  compliance  with  specific  performance  and  serviceability  criteria. 
If  the  component  has  failed  during  that  time  period,  its  life  is  over.  The  model  then 
computes  the  component’s  lifetime  and  dixrability  costs.  If  the  structure  has  not  failed, 
but  instead  a  repair  has  been  triggered,  then  the  repair  is  executed,  the  durability  cost 
revised,  and  the  element  is  ready  to  enter  the  next  time  period.  If  the  element  meets 
all  criteria  then  it  simply  passes  to  the  next  time  period.  The  simulation  continues 
until  the  design  life  has  expired  or  the  component  has  failed. 

Clearly,  the  success  of  this  durability  model  depends  upon  how  well  the  environmental 
factors  can  be  simulated  and  how  well  the  damage  processes  can  be  modeled  and 
repair  strategies  defined.  However,  the  paradigm  can  be  implemented  with  rather 
primitive  environmental  and  damage  models  providing  that  more  sophisticated  models 
can  be  put  in  place  of  the  simple  ones  as  they  become  available.  In  this  work  some 
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Figures  are  found  at  the  end  of  their  corresponding  chapters. 
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fairly  simple  models  of  damage  are  adopted  to  demonstrate  the  S5mergy  of  the  coupled 
approach. 


Environmental  Stimuli 

Over  its  lifetime,  a  material  component  will  be  subjected  to  many  environmental 
factors,  and  these  factors  will  generally  fluctuate  with  time.  In  this  model,  the  BMDM 
assume  that  the  changes  in  these  environmental  stimuli  over  time  are  known.  Three 
environmental  effects  will  be  considered  here:  the  ambient  temperature  T(0,  the 
ambient  moisture  6(0  due  to  humidity  and  precipitation  (expressed  as  a  fraction  of 
saturation),  and  the  concentration  of  aggressive  chemicals  or  agents  of  decay  c{t) 
(expressed  as  a  fraction  of  the  maximum  concentration).  Different  environments  can 
be  modeled  by  assigning  different  histories  for  these  three  stimuli.  Two  h5T)othetical 
environments  will  be  considered:  tropical  (possibly  indicative  of  Guam)  and  conti¬ 
nental  United  States  (possibly  indicative  of  Chicago). 

Temperature,  humidity,  and  chemical  concentration  generally  fluctuate  according  to 
daily  and  seasonal  changes.  For  the  current  purpose,  the  monthly  averages  of  these 
stimuli  are  used.  When  daily  fluctuations  are  important,  as  they  would  be  in  a  freeze- 
thaw  model,  they  will  be  accounted  for  explicitly  in  damage  models.  Regular  seasonal 
changes  will  be  augmented  with  random  fluctuations  to  simulate  the  variability  of 
climatic  effects. 

The  fluctuation  of  temperature  is  given  by  the  expression 

T(t)  =  To+ATsin(nt/6)  +6Tn(t)  [Eq  1] 

where  T^  is  the  average  yearly  temperature,  AT  is  half  the  difference  between  the 
highest  average  monthly  temperature  and  the  lowest  average  monthly  temperature, 
6T  is  the  average  fluctuation  in  temperature  in  1  month,  and  T|(i)  is  a  uniformly 
distributed  random  number  between  -1  and  1.  The  time  t  is  measured  in  months.  The 
continental  U.S.  temperature  model  has  Tq  =  40°F,  AT  =  50°F,  and  6T  =  10“F.  The 
tropical  model  has  Tq  =  85“F,  AT  =  10®F,  and  6T  =  5“F. 

The  fluctuation  in  moisture  is  given  by  a  similar  expression, 

0(t)  =  eo+Aesin(nt/6)+6en(t)  [Eq2] 

where  %  is  the  average  annual  moisture  content,  A0  is  half  the  difference  between  the 
highest  average  monthly  moisture  content  and  the  lowest  average  monthly  moisture. 
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60  is  the  average  fluctuation  in  moisture  in  1  month,  and  iiCO  is  a  uniformly 
distributed  random  number  between  -1  and  1.  The  continental  U.S.  moisture  model 
has  00=  0.5,  A0  =  0.3,  and  60  =  0.05.  The  tropical  model  has  0o=  0.9,  A0  =  0.05,  and 
60  =  0.01. 


The  chemical  environment  also  can  vary  by  season,  but  for  current  purposes  it  is 
modeled  simply  as  a  constant  average  value  with  random  fluctuations, 

c(t)  =  c,+6cn(t)  [Eq3] 

where  Cq  is  the  average  concentration  of  aggressive  chemicals,  6c  is  the  expected 
fluctuation  in  that  average,  and  r|(0  is  a  uniformly  distributed  random  number 
between  -1  and  1.  The  continental  U.S.  model  has  Cq=  0.7  and  6c=  0.3.  The  tropical 
model  has  Co=  0.2  and  6c=  0.1. 


General  Framework  for  Damage  Models 

Damage  mechanisms  for  various  materials  are  generally  determined  by  factors  that 
occur  at  a  microscopic  scale.  Many  of  these  processes  are  still  poorly  understood.  It 
is  beyond  the  scope  of  the  current  investigation  to  attempt  to  model  damage  at  the 
microscale.  One  might  eventually  replace  each  of  the  various  damage  models  used 
here  with  more  accurate  ones,  but  for  current  purposes  simple  models  are  used  to 
qualitatively  reflect  the  physical  processes  that  are  taking  place. 

Let  p(t)  represent  a  t3q)ical  damage  parameter.  Depending  upon  the  mechanism  of 
damage,  the  evolution  of  this  damage  parameter  will  be  governed  by  one  of  three  basic 
damage  models:  linear  growth,  exponential  growth,  or  polynomial  growth. 

Linear  Growth  Model 

A  linear  growth  model  does  not  depend  on  the  current  state  of  damage,  but  has  a  rate 
of  growth  governed  by  some  function  of  the  environment: 


dt 


=fCr.0.c) 


[Eq4] 


P(0)=Po 
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For  this  model,  if  the  environment  is  constant,  damage  grows  linearly  with  time: 

p(t)  =  Po*<(T,8.c)t  [Eq5] 

Scaling  of  concrete  caused  by  freeze-thaw  cycling  is  an  example  of  linear  growth 
damage. 

Exponential  Growth  Model 

In  an  exponential  growth  model,  the  rate  of  change  in  the  damage  parameter  depends 
directly  upon  the  amount  of  damage  currently  present  (sometimes  raised  to  some 
power  m): 


5=fCr.e.c)p(t) 

at 


P(0)=Po 


[Eq6] 


For  this  model,  if  the  environment  is  constant,  the  damage  evolves  exponentially  with 
time  as  follows; 


p(t)  = 


[Eq7] 


The  solution  to  the  exponential  growth  model  for  Po=l  and  constant  forcing  function 
f=B  is  shown  in  Figure  2  for  different  values  of  B.  Damage  to  the  surface  protection 
provided  by  paint  is  modeled  as  an  exponential  growth  process. 

Polynomial  Growth  Model 

A  polynomial  growth  model  represents  a  diffusion  process.  The  rate  of  change  of  the 
damage  parameter  is  inversely  proportional  to  the  amount  of  damage  present  (often 
raised  to  some  power  m): 


dp  _  f(T,9.c) 
dt  '  p(t) 


[EqS] 


P(0)=Po 
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For  this  model,  if  the  environment  is  constant,  the  damage  evolves  as  a  pol)momial 
with  time  as  follows: 

p(t)  =  /po^2f(T.e.c)t  IEq9] 

One  of  the  most  important  features  of  this  model  is  that  it  is  usually  associated  with 
failure  in  a  finite  amount  of  time.  The  solution  to  the  exponential  growth  model  for 
Po=l  and  constant  forcing  fimction  f^B  is  shown  in  Figure  3  for  different  values  of  B. 
Corrosion  of  steel,  decay  of  wood,  and  static  fatigue  of  concrete  are  examples  of  the 
diffusion  (polynomial  growth)  model. 

In  each  case,  the  environmental  function  driving  the  process  will  be  different.  In  most 
cases  a  four-parameter  model  was  selected,  with  one  tunable  parameter  for  each  of  the 
three  environmental  effects,  and  one  to  accurately  represent  the  overall  rate  of 
damage.  Some  of  these  functional  dependencies  come  from  first  principles,  e.g., 
Arrhenius’s  Law  (Arrhenius  1907)  is  used  to  characterize  dependence  upon  tempera¬ 
ture;  others,  however,  are  simply  plausible  heuristic  models,  e.g.,  the  exponential 
models  often  used  for  moisture  and  chemical  concentration.  A  typical  example  of  one 
of  these  functions  is  given  by: 

f  (T,e.c)=D,te(t)]"[1  +c(t)]"’e  -P/T(t)  [Eq  1 0] 


This  type  of  model  “does  the  right  thing”  in  the  sense  that  higher  moisture  content 
drives  damage  faster,  but  there  is  a  premium  on  moisture  contents  near  saturation, 
0=1,  and  damage  stops  if  the  element  dries  out.  (Not  all  models  have  this  feature.) 
Higher  concentrations  of  aggressive  chemicals  speed  up  the  rate  of  damage,  but 
damage  will  occur  even  if  no  chemicals  are  present.  High  temperatures  speed  up 
damage;  low  temperatures  nearly  stop  the  damage  process. 

Numerical  Integration  of  the  Damage  Models 

The  fimctions  /(T,  6,  c)  driving  the  damage  processes  will,  in  general,  not  be  constant 
with  time.  Thus,  the  rate  equations  must  be  integrated  numerically.  For  current 
purposes  the  generalized  midpoint  rule  is  used.  The  incremental  linear  growth  model 
of  Equation  4  then  takes  the  form 
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P(Vi)=P(t|)+f(To,.ee„CjAt 


P(0)=Po 


[Eqll] 


where  and  the  environmental  factors  are  a  weighted  average  of  the  values 


at  and  as  follows: 


c„  =  (1-a)c(t|)+ac(th,i) 


0„=(1-a)e{t^+ae(t|,i) 


[Eq12] 


T„H(1-a)T(t,)+aT(t,.i)  ] 

In  general,  one  should  use  a=0.5,  i.e.,  the  trapezoidal  rule.  The  other  equations  are 
integrated  similarly.  The  exponential  growth  model  has  the  incremental  form 


P(t|.i)=P(t|)e 


P(0)=Po 


f(T„,e,.<s„)At 


[Eq  13] 


The  polynomial  growth  model  has  the  incremental  form 


p(V,)=«h2f(TX^ 


P(0)=Po 


[Eq  14] 


Throughout  this  work  details  of  computation  are  given  only  for  exceptions  to  the  above 
cases  and  for  computations  that  help  determine  the  coefficients  of  the  models. 


Cost  Modeling 


The  user  can  specify  an  MER  strategy  in  accord  with  serviceability  requirements, 
safety  factors,  or  other  factors  edfecting  the  performance  of  each  component.  Each 


16 


USACERL  TR  96/24 


material  will  be  subject  to  different  constraints,  but  for  current  purposes  the  strategy 
is  modeled  as  follows:  the  vth  repair  act  costs  C^,  dollars  and  restores  the  element 
properties,  as  well  as  certain  of  the  model  parameters,  to  a  fraction  of  the  most 
recent  repaired  value  or  to  a  fraction  of  its  original  value.  An  element  that  fails 
during  a  time  step  must  be  replaced — ^it  cannot  be  repaired. 

The  quality  of  a  component  or  a  repair  job  will  depend  on  many  factors.  For  example, 
the  quality  of  a  paint  job  depends  upon  the  quality  of  the  paint,  the  quality  of  surface 
preparation,  and  the  quality  of  paint  application.  To  simplify  this  discussion  it  is 
hypothesized  that  the  overall  quality  can  be  indexed  with  a  single  value  based  on  an 
assumption  that  a  person  who  buys  the  best  paint  would  also  be  inclined  to  do  a  good 
job  preparing  the  surface  and  applsdng  the  paint.  In  other  words,  this  model  simply 
cannot  account  for  a  person  doing  a  sloppy  job  with  good  paint.  The  quality  q  will  be 
specified  as  a  number  between  0  (poorest  quality)  and  1  (best  quality). 

The  cost  associated  with  making  a  repair  is  also  a  function  of  the  quality  of  that  repair. 
Clearly,  the  best  quality  repair  will  cost  the  most  and  the  worst  quality  repair  the 
least.  The  cost  of  repair  is  modeled  as 


C(q)=— ^  [Eq15] 

l-q(l-p) 

where  Cq  is  the  cost  of  the  poorest  quality  repair  and  p  is  the  ratio  of  the  cost  of  the 
poorest  quality  repair  to  the  cost  of  the  best  quality  repair.  The  model  3delds  a  “pro¬ 
gressive  tax”  on  quality  as  it  nears  perfection,  but  is  relatively  insensitive  to  cost  at 
lesser  qualities. 

The  user  must  specify  the  initial  cost  of  the  component,  to  facilitate  comparison  among 
materials  as  well  as  the  cost  of  replacement.  In  some  models  certain  things  can  be 
repaired  while  others  cannot.  In  those  cases,  the  repair  strategy  applies  only  to  items 
that  can  be  repaired.  The  maintenance  scheme  requires  the  specification  of  the  least- 
expensive  repair,  the  most  expensive  repair,  the  initial  quality  qo  of  the  component, 
and  the  quality  of  any  subsequent  repair  q„  (as  well  as  the  criteria  that  determine  the 
time  of  repair).  These  specifications  determine  the  durability  cost. 

The  time  value  of  money  is  considered  by  finding  the  present  value  of  a  future  amoimt 
according  to  the  following  formula  for  single  payment  present  worth  (SPPW): 


SPPW=(1+i)-" 


[Eq  16] 
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where  i  is  the  effective  discount  rate  per  period  (1  month  for  the  current  simulations) 
and  n  is  the  number  of  periods.  If  the  cost  of  a  repair  is  at  the  time  of  repair,  then 
the  value  of  that  repair  at  time  zero  is 

C,=C,(l+i)-"  [Eq17] 

where  n  is  the  number  of  periods  elapsed  since  time  zero.  For  the  accumulation  of 
durability  costs,  all  costs  are  referred  back  to  time  zero.  Therefore,  a  material 
requiring  repair  in  40  years  (for  example)  can  be  considered  less  costly  in  terms  of 
current  dollars  than  a  material  requiring  repair  in  30  years.  The  time  value  of 
money — ^that  is,  the  idea  that  a  dollar  today  is  inherently  more  valuable  than  the  same 
dollar  next  year — ^may  play  an  important  role  in  comparing  a  high-maintenance 
material  (such  as  steel)  with  a  low-maintenance — ^but  initially  more  expensive — 
material  (such  as  aluminum). 


18 


USACERL  TR  96/24 


Figure  1.  Durability  model  concept. 
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Figure  2.  Exponential  growth  equations  with  constant  coefficients. 
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3  Modeling  the  Durability  of  a  Structural 
Element 

Description  of  Model  Problems 

As  a  model  problem,  consider  the  durability  of  a  beam  subjected  to  a  loading  Pit)  at  its 
midspan,  as  shown  in  Figure  4.  The  materials  available  for  beam  construction  are 
reinforced  concrete,  steel,  wood,  aluminum,  or  fiber-reinforced  plastic  (FRP).  Over  its 
lifetime,  the  beam  is  subjected  to  environmental  factors  that  fluctuate  with  time.  As 
noted  in  Chapter  2,  this  problem  focuses  on  three  factors: 

•  ambient  temperature  T(f) 

•  ambient  moisture  6(f)  due  to  hiimidity  and  precipitation,  expressed  as  a  fraction 
of  saturation 

•  concentration  of  aggressive  chemicals  or  agents  of  decay  c(f),  expressed  as  a 
fraction  of  the  maximum  concentration. 

For  the  simulations  P(t)  takes  the  following  form: 

P(t)=P„+6Pn"(t)  [Eq18] 

where  Pq  is  the  expected  load  and  bP  is  the  amplitude  of  the  fluctuation  in  the  load. 
As  before,  ri(<)  is  a  uniformly  distributed  variate  (random  variable). 

The  beam  must  continue  to  function  according  to  its  intended  purpose,  that  is,  carry 
the  loading  Pit)  without  excessive  deflection.  For  each  material  there  are  strength  and 
serviceability  requirements  in  terms  of  cross-sectional  properties  hke  the  area  A(t)  and 
the  moment  of  inertia  I(t),  which  are  degrading  through  various  damage  mechanisms. 
Both  the  onset  of  failure  (or  dysfunction)  and  the  M&R  cost  are  considered. 
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Strength  Design  Criterion 

It  is  assumed  that  the  beam  is  designed  to  have  adequate  bending  strength  at  time 
^=0.  The  strength  of  the  beam,  as  a  function  of  time,  is  designated  by  Pf^t).  The  initial 
strength  of  the  beam  is  thus  PuiO)  and  meets  the  design  criterion 

Pu(0)"<l>oPe.^  [Eq19] 

where  (|)o  is  the  initial  factor  of  safety  and  P^p,,ted  is  the  maximum  expected  applied 
load.  As  the  properties  of  the  beam  degrade  over  time  the  strength  will  decrease, 
thereby  compromising  the  factor  of  safety.  The  current  design  factor  of  safety  is  given 
by 


Pu(t) 

[Eq20] 

and  define  the  current  actual  factor  of  safety  is  given  by 

P  (t) 

IEq21] 

The  current  design  factor  of  safety  will  be  used  as  an  index  for  repair  decisions  (e.g., 
repair  if  4)^  <  0.8(|)o)  while  the  current  actual  factor  of  safety  will  be  used  to  detect 
failure  of  the  system,  i.e.,  (|)a  <  1. 

Serviceability  Design  Criterion 

The  main  serviceability  design  criterion  is  the  guard  against  excessive  deflection.  The 
initial  design  will  have  deflection  limits  for  a  variety  of  reasons,  many  of  which  will 
have  little  to  do  with  the  long-term  performance  of  the  beam.  As  a  consequence,  a 
long-term  limit  to  midspan  deflection  is  simply  specified  as  Wi^.  The  current 
serviceability  factor  of  safety  is  given  as 


w(t) 


w» 


[Eq22] 


The  current  serviceability  factor  of  safety  will  be  used  as  an  index  for  repair  decisions 
(e.g.,  repair  if  (J),  >  1)  or  to  detect  functional  failure  of  the  system  (e.g.,  the  component 
is  dysfunctional  if  (l)^  >  10). 
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Durability  Model  for  a  Reinforced  Concrete  Beam 


Degradation  of  the  reinforced  concrete  beam  is  caused  by  corrosion  of  the  reinforce¬ 
ment  bars,  static  fatigue  of  the  compressive  strength,  and  scaling  of  the  top  surface, 
as  shown  schematically  in  Figure  5.  The  strength  of  the  reinforced  concrete  beam  at 
time  t  is  given  by 


h(t)-0.59 


A.(t)fy 

Oo(t)b 


[Eq23] 


where  A/^)  =  n^izlrQ  -  is  the  area  of  steel  bars  corroded  to  a  depth  ofac(/),  h(t)  = 
ho  -  sit)  is  the  beam  depth  scaled  to  a  depth  sit),  and  aj,t)  is  the  concrete  compressive 
strength.  The  3deld  strength  of  the  steel,  fy ,  and  the  width  of  the  beam,  h,  are  taken 
as  constants. 

Although  reinforced  concrete  is  subject  to  long-term  creep  deflections,  they  are  not 
considered  in  this  problem. 

Corrosion  of  Reinforcement  Bars 

The  beam  is  reinforced  with  steel  bars  of  initial  radius  r^.  Due  to  the  history  of 
temperature  T,  moisture  0,  and  corroding  chemicals  c,  the  bars  have  a  current  radius 
of  r,  reduced  from  the  original  radius  by  the  amount  x.  Figure  6  shows  the  measure 
of  material  loss  caused  by  corrosion. 

Corrosion  of  the  reinforcement  bars  will  depend  on  how  much  moisture  or  other 
chemical  agents  are  present.  In  lieu  of  a  more  sophisticated  model  of  moisture 
transport,  it  is  assumed  that  the  concrete  has  pre-existing  cracks  that  allow 
environmental  moisture  unimpeded  access  to  the  reinforcement  bars.  The  following 
rate  equation  describes  the  loss  of  bar  radius  caused  by  corrosion: 


dx  ^ 

dt  x(t) 

x(0)  =  0 

The  driving  function  takes  the  form 


[Eq24] 


f,(r,e,c)  =  D„[9(t)ni+c(t)re-'>/"® 


[Eq25] 
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where  0(^)  is  the  environmental  moisture  content  as  a  percentage  of  saturation,  c(t)  is 
the  concentration  of  the  chemical  environment,  and  T(t)  is  the  temperature  (in  degrees 
Kelvin).  The  parameters  n,  m,  Dq,  and  p  are  material  constants  that  will  be  fit  to 
anecdotal  data;  in  other  words,  all  other  factors  will  be  held  constant  while  the 
published  variation  caused  by  a  single  effect  are  fit. 

The  solution  to  the  rate  can  be  obtained  by  numerical  quadrature  (squaring)  to  give 
the  loss  of  radius  at  time  in  terms  of  the  loss  at  time  augmented  by  the  corrosion 
that  takes  place  between  those  two  times: 


x(V,)  =  \/x"(^)-2Atf,(T„,e„.c„) 


x(0)=0 


[Eq26] 


The  main  idea  behind  this  corrosion  model  is  as  follows:  because  corrosion  is  mainly 
the  result  of  exposure  to  moisture,  the  rate  of  corrosion  should  increase  as  the  moisture 
content  increases — and  should  effectively  stop  in  the  absence  of  moisture.  Further¬ 
more,  nearly  dry  conditions  should  be  considerably  less  corrosive  than  nearly  satu¬ 
rated  conditions.  Therefore,  in  Equation  25,  the  exponent  n  is  greater  than  1.  The 
greater  the  exponent,  the  more  accentuated  this  effect  will  be. 


Corrosion  is  exacerbated  by  the  presence  of  aggressive  chemicals,  but  it  will  occur  even 
in  the  absence  of  chemicals.  The  greater  the  concentration  of  aggressive  chemicals,  the 
faster  the  corrosion.  Therefore,  the  exponent  m>l.  The  greater  the  exponent,  the 
more  accentuated  this  effect  will  be. 


According  to  Arrhenius’s  Law  (Arrhenius  1907),  the  rate  of  corrosion  depends  upon  the 
temperature.  The  material  will  corrode  faster  at  a  higher  temperature.  The  coefficient 
P  controls  the  rate  of  corrosion  caused  by  an  increase  in  temperature.  Corrosion  would 
be  expected  to  drop  to  arbitrarily  low  rates  as  the  temperature  drops  below  the 
freezing  point,  and  would  occur  fastest  at  the  top  of  the  temperature  range  being 
considered.  The  maximum  enAdronmental  temperature  may  be  designated  as  T^.  If, 
at  any  given  instant  the  rate  of  corrosion  at  T^  were  y  times  the  rate  at  eT^  then  the 
exponent  p  must  be  given  by 


3  =  -^T.InY 


[Eq27] 
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For  example,  if  the  rate  of  corrosion  is  one-fourth  as  fast  at  half  the  maximum 
temperature  of  325  °K,  thenp  =  450  °K. 

Static  Fatigue  In  Concrete 

At  loads  that  cause  the  compressive  stress  to  approach  the  strength  of  concrete, 
damage  is  incurred  through  the  mechanism  of  static  fatigue.  This  phenomenon  results 
in  a  decrease  in  the  strength  of  concrete  over  time.  Let  Oo(^)  be  the  concrete  strength 
and  o(t)  (positive  compressive)  the  applied  stress  at  time  t.  Assuming  that  the 
material  exhibits  elastic  behavior,  one  can  relate  the  stress  to  the  applied  load  though 
the  relationship  a{t)  —  Pit)/S,  where  S  is  the  section  modulus  of  the  beam.  For  current 
purposes  it  is  hypothesized  that  the  strength  degrades  according  to  a  power  law  as 
follows 


o(t) 

Oo(t)-o(t) 


[Eq28] 


Oo(0)=O,  j 

where  B  and  n  are  the  parameters  describing  the  material.  Assuming  that  the  applied 
stress  remains  constant  from  time  to  ^i+l>  Equation  28  can  be  integrated  to  give  the 
strength  at  the  current  time  as 


=  o(t,)+  [[Oo(V-o(t,)l-^  -(n+1)Bo,o"(VAt]  IE*!  29] 

This  model  will  predict  a  decrease  in  strength  with  time,  but  it  can  be  rather  sensitive. 
Therefore,  one  must  take  care  to  check  for  material  failure  within  a  time  step.  The 
time  to  failure,  expressed  as  a  fraction  of  the  current  time  increment,  is  given  as 


[Oo(t,)-o(t,)r^^ 

(n+1)Bd„o"(t,)At 


[Eq30] 


If  !/•  <  1,  then  the  material  has  failed  during  the  current  time  increment.  Because  one 
must  compute  time  to  failure  at  every  step  an3rway.  Equation  29  can  be  simplified  to 
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1 


[Eq31] 


oJP)  =  0„ 


The  constants  can  be  fit  to  data  obtained  from  tests  in  which  the  time  to  failure  is 
recorded  for  constant  stress  values.  In  such  an  instance,  failure  is  defined  as  Opitf)  = 
o,  and  time  to  failure  is  given  by 


(n+1)Bo<,o" 


[Eq32] 


Assume  that  two  tests  are  available.  Let  (o^ ,  )  be  the  applied  stress  and  time  to 

failure  of  the  first  test,  and  let  (a^,  be  the  applied  stress  and  time  to  failure  of  the 
second  test.  The  exponent  of  the  model  can  be  determined  from  these  two  tests  by  the 
following  relationship; 


nln 


Ol(Oo-02) 


OaK-Oi) 


=  ln 


<2(Po-Oi) 

ti(Oo-02) 


[Eq33] 


After  n  is  determined,  B  can  be  determined  from 


(n+1)d„o?ti 

or 


[Eq34] 


(n+1)o„<^ 

For  example,  if  at  =  0.90^  with  ti  =  0.484  and  O2  =  0.75oowith  ^2  =  294  then  one 
obtains  n  =  5  and  JB  =  3.5  x  10  ®. 


Concrete  Seating 

Scaling  of  a  concrete  surface  is  caused  largely  by  freeze-thaw  cycling  of  the  environ¬ 
ment.  The  mechanism  causing  material  damage  from  freeze-thaw  cycling  is  extremely 
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complex,  involving  the  moisture  content  at  the  time  of  freezing,  the  rate  of  freezing, 
and  the  pore  structure  of  the  material  (which  cjm  also  alter  the  freezing  temperatures). 
To  develop  a  simplified  model  that  proves  the  concept  by  making  qualitatively  correct 
projections — ^which  is  the  present  purpose — ^it  is  necessary  to  make  some  broad 
assumptions. 

First  it  is  assumed  that  the  number  of  freeze-thaw  cycles  is  related  to  the  average 
temperature  multiplied  by  the  time  over  which  the  temperature  was  averaged. 
(Clearly  the  model  will  behave  best  if  one  takes  time  increments  small  enough  to  sense 
seasonal  changes,  but  with  this  model  one  does  not  need  to  take  them  small  enough 
to  resolve  daily  fluctuations).  Next,  it  is  assumed  that  the  surface  deteriorates  from 
the  outside-in,  with  the  amount  of  deterioration  during  a  time  step  being  given  by  the 
number  of  freeze-thaw  cycles.  The  depth  of  scaling  s(0j  shown  in  Figure  7,  is  thus 
modeled  as 


^=f,(T.0.c) 

s(0)=0 

where  the  driving  function  has  the  form 

f3(r.e,c)  =  C<T„-T(t)>[0(t)r[1  +c(t)r 


[Eq351 


[Eq36] 


and  where  the  parameters  C,  n,  and  m  affect  the  rate  of  scaling  caused  by  moisture 
and  concentration  of  chemicals,  much  like  the  corrosion  model.  The  notation  <x> 
means  that  <jc>=  x  if  x>0  and  <x>=  0  if  x<0.  The  value  Tq  is  the  fixed  temperature 
above  which  no  freeze-thaw  cycling  is  expected  to  occur.  For  example,  if  the  average 
temperature  over  a  month  is  10  °F  one  might  expect  daily  fluctuations  to  result  in 
freeze-thaw  cycling.  While  one  might  expect  less  cycling  for  very  cold  temperatures, 
it  is  assumed  for  the  current  problem  that  the  temperature  of  the  material  is  affected 
by  its  ability  to  absorb  the  sun’s  heat.  For  example,  blacktop  can  thaw  even  if  the  air 
temperature  does  not  exceed  the  freezing  temperature  for  water. 

Assiiming  the  environmental  factors  to  be  constant  over  the  time  step.  Equation  35  can 
be  integrated  to  give 
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S(U=S(V^A«en'a.6a.C„) 

S(0)=0 


[Eq37] 


The  model  parameters  can  be  adjusted  to  fit  anecdotal  data.  The  influence  of  moisture 
should  decrease  rapidly  as  the  moisture  level  falls  below  satiiration.  Thus,  an 
exponent  like  n=5  might  be  reasonable.  The  influence  of  the  chemical  concentration 
might  be  adequately  reflected  by  a  value  like  m=3.  The  value  of  C  might  best  be 
determined  by  considering  the  time  it  would  take  to  scale  to  a  certain  depth  in  the 
harshest  environment. 


Durability  Model  for  Steel 

Degradation  of  the  steel  beam  is  caused  by  corrosion  of  the  exposed  surfaces,  as  shown 
schematically  in  Figure  8.  The  strength  of  the  steel  beam  at  time  t  is  given  by 

Pu(t)  =  ^[h%-2x(t)]  +  4hb[t,-2x(t)]]  [Eq  38] 

where  all  of  the  plate  faces  are  corroded  to  a  depth  of  x(t).  The  strength  of  the  steel  fy, 
the  depth  of  the  beam  h,  and  the  width  of  the  beam  b  are  taken  as  constants.  It  is 
assumed  that  there  are  no  long-term  deflection  problems  for  a  steel  beam. 

Corrosion  of  Steel 

The  steel  beam  has  an  I-t3q)e  cross-section  with  depth  h  and  flange  width  b.  It  has  a 
flange  thickness  tf  and  a  web  thickness  that,  due  to  the  history  of  temperature  T, 
moisture  0,  and  corroding  chemicals  c,  have  corroded  by  amoimt  x  on  each  face.  Figure 
9  shows  the  measure  of  material  loss  caused  by  corrosion. 

Corrosion  of  the  steel  beam  will  be  modeled  in  the  same  manner  as  the  reinforcement 
bars  were  modeled  for  the  reinforced  concrete  beam,  but  a  modification  must  be  made 
in  the  environmental  exposure.  Since  steel  is  almost  always  painted,  a  paint  model 
is  considered  in  conjunction  with  the  corrosion  model.  The  moisture  and  concentration 
of  aggressive  chemicals  in  the  corrosion  model  will  simply  be  the  amount  transported 
through  the  surfacing  material. 
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The  following  diffusion  equation  describes  the  loss  of  plate  thickness  caused  by 
corrosion  after  the  vth  repair: 


x°(0)  =  0 


dx^  fcfT.ep.Cp) 
dt  x''(t) 

x''(0=x''-’(t:-') 


[Eq39] 


where  the  vth  repair  is  executed  at  time  to"  s  with  the  initial  condition  for  the  next 
stage  being  given  by  the  end  condition  of  the  previous  stage.  The  cross-section  starts 
in  an  rmdamaged  condition.  The  corrosion-driving  function  is  given  by  the  expression 


f.(T,ep,Cp)  =  Djep(t)r[UCp(t)re 


[Eq40] 


where  0p«)  is  the  moisture  content  as  a  percentage  of  saturation,  and  Cpit)  is  the 
concentration  of  aggressive  chemicals  transported  through  the  surfacing  material.  T(t) 
is  the  absolute  temperature  (Kelvin).  The  parameters  n,  m,  D„,  and  p  govern  the  rate 
of  corrosion.  These  parameters  can  be  established  in  exactly  the  same  manner  as  in 
the  corrosion  model  for  reinforcement  bars  in  the  concrete  beam  model. 

Surface  Porosity  Model 

The  porosity  of  the  surface  pit)  is  governed  by  a  damage  model  that  addresses  both 
coating  and  substrate  condition.  The  rate  of  change  in  porosity  is  assumed  to  be 
proportional  to  the  amount  of  porosity,  and  is  driven  by  a  function  of  the  environmen¬ 
tal  factors  as  follows: 


^  =  -fp(T,e,c,v)p''(t) 


P''(0  =  9p(qv) 


[Eq41] 


p»(0)  =  1 
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where  the  function  driving  the  change  in  porosity  is 

fp(T.0,c.v)  =  Dp[T(t)f[1  +0{t)f  [1  +c(t)r  [Eq  42] 


and  where 


T(t)  = 


t(t)-Tp? 

Ti 


tEq431 


is  a  normalized  measure  of  the  temperatxire.  It  has  been  noted  that  damage  to  paint 
generally  occurs  at  a  greater  rate  for  very  cold  temperatures  and  for  very  hot 
temperatures,  but  slows  for  intermediate  temperatures  (Tooke  1980).  For  this  model 
the  reference  temperatures  T„  =  300  ®K  and  Tj  =  20  “K  are  used.  Deterioration  speeds 
up  in  the  presence  of  moisture  and  aggressive  chemicals,  but  it  will  also  proceed  in  a 
dry,  chemically  neutral  environment. 

The  quality  of  the  surfacing  determines  the  condition  to  which  the  porosity  index 
returns.  The  function  used  in  this  model  is  simply  taken  to  be 

gp(qv)  =  (qv)® 

The  parameters  n,  m,  and  k  are  taken  to  be  functions  of  the  quality  of  the  surfacing 
(i.e.,  material  used,  surface  preparation,  etc.).  The  user  must  specify  the  quality  q  of 
the  initial  surfacing  and  of  any  subsequent  repair  surfacing.  The  quality  q  will  be 
specified  as  a  number  between  0  (poorest  quality)  and  1  (best  quality).  Then  the 
parameters  n  and  m,  at  the  vth  repeur,  can  be  specified  in  accord  with  the  quality  of 
the  repair  using  the  following  relationships 


n„  =  4/q 


iTV  =  3/q 


[Eq45] 


K,  =  1/q 


The  maintenance  strategy  then  consists  of  specifying  the  times  Bind  qualities  of  the 
resurfacing. 
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Equation  41  can  be  numerically  integrated  to  give  the  following  value  of  the  porosity 
at  time  in  terms  of  the  porosity  at  time  after  the  vth  repair: 


p''(0=gp(qv) 


[Eq46] 


p‘>(0)  =  1  ] 

where  A<  =  and  the  environmental  factors,  as  before,  are  the  trapezoidal 

averages  of  values  at  t''i  and 

As  the  porosity  index  increases  with  damage,  the  resistance  of  the  surfacing  material 
to  moisture  and  chemical  transport  decreases.  The  following  simple  transport  model 
gives  the  moisture  and  chemical  concentration  at  the  surface  of  the  steel  in  terms  of 
the  environmental  moisture  and  chemical  concentration: 


epW^i-PWlt) 


[Eq47] 


Cp(t)  =  [1-P''(t)]c(t) 


The  solution  to  the  rate  equation  for  corrosion  can  now  be  obtained  by  numerical 
quadrature  to  give  the  loss  of  plate  thickness  at  time  in  terms  of  the  loss  at  time 

augmented  by  the  corrosion  that  takes  place  between  those  two  times: 


x(Ci)  =  v'x2(t,'')+2Atf,(T„,ep„.Cp„) 


x“(0)  =  0 


[Eq48] 


x''(0=x''-’(O 


where  0p„  and  c^^are  evaluated  in  a  manner  analogous  to  0p  and  Cp  in  Equation  12. 
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Durability  Model  for  Wood 

Degradation  of  the  wood  beam  is  caused  by  decay  of  the  exposed  surfaces,  as  shown 
schematically  in  Figure  10.  The  strength  of  the  wood  beam  at  time  t  is  given  by 


Pu(t)  =  ^[h-2x(t)l2[b-2x(t)]  [Eq49] 

where  all  of  the  beam  faces  are  decayed  to  a  depth  of  a:(^).  The  modulus  of  rupture  of 
the  wood  the  depth  of  the  beam  h,  and  the  width  of  the  beam  b  are  constants. 

Decay  of  Wood 

The  wood  beam  has  a  rectangular  cross-section  with  depth  h  and  width  b  that,  due  to 
the  history  of  temperature  T,  moisture  0,  and  corroding  chemicals  c,  have  decayed  by 
the  amount  x  on  each  face.  Figure  11  shows  the  measure  of  material  loss  caused  by 
rotting. 

The  following  diffusion  equation  describes  the  loss  of  material  caused  by  decay: 


dx  _  tdO’iQ) 
dt  x(t) 

x(0)  =  0 


[EqSO] 


The  cross-section  starts  in  an  undamaged  condition.  The  decay-driving  function  is 
given  by  the  expression 


fd(f.e)  =  Cd[1-0-^“>][0(t)r 


[Eq51] 


where 


T(t). 


T(t)-TJ^ 

Ti 


[Eq52] 


is  a  normalized  measure  of  the  temperature.  It  has  been  noted  that  decay  generally 
occurs  only  over  a  certain  range  of  temperatures  (American  Institute  of  Timber 
Construction  1974).  For  this  model  the  reference  temperatures  are  taken  as  =  300 
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°K  and  Tj  =  20  °K.  The  deterioration  speeds  up  in  the  presence  of  moisture  but  will 
not  take  place  in  a  dry  environment. 

The  parameters  n,  m,  and  k  depend  on  the  quality  (hardness)  of  the  wood  Q,  with  Q=1 
being  the  softest  wood  and  Q=10  being  the  hardest  wood.  The  following  values  are 
recommended: 


n  =  ,/Q 


[Eq53] 


As  before,  the  solution  to  the  rate  can  be  obtained  by  numerical  quadrature  to  give  the 
loss  of  material  at  time  h+I  in  terms  of  the  loss  at  time  ti  augmented  by  the  decay  that 
takes  place  between  those  two  times: 


[Eq54] 


x(0)=0  J 

where  At  =  -  ti  and  the  environmental  functions,  as  before,  are  the  trapezoidal 

average  of  values  at  and  t^^j. 


Creep  Deflection  of  Wood 

Wood  is  subject  to  creep  deflections  caused  by  long-term  loading.  To  model  the 
long-term  creep  behavior  a  Norton  exponential  creep  law  (Kachinov  1986)  is  employed. 
The  law  is  assumed  to  apply  to  the  moment-curvature  relationship  of  the  beam  as 
follows: 


ii(z,t)  =  ^+BM'"(z,t)  [Eq55] 

where  k  is  the  rate  of  change  of  curvature  of  the  beam,  M  is  the  bending  moment,  and 
El  is  the  flexural  modulus,  which  may  vary  with  time  due  to  loss  of  cross-section.  The 
parameters  B  and  m  are  the  creep  coefficients. 
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For  the  model  problem  M(z)  -  Pz  1 2,  enabling  the  integration  of  Equation  55  with 
respect  to  the  spatial  dimension  z,  and  resulting  in  an  equation  for  the  rate  of  change 
of  deflection  as  follows; 


w(z,t)  =  (KZitlPlb+vlzlP"*!!)  [Eq  56] 

where  the  beam  deflected-shape  functions  are  given  by 


<l>(z,t)  = 


3zL^-4z^ 

48EI(t) 


^  p(mf2)zL'"-^-2'"-^z'"-^ 
(m+1)(m+2)22"*’ 


[Eq57] 


The  maximum  deflection  occurs  at  z=LI2.  Any  variable  denoted  with  a  hat  ('^)  both 
depends  on  z  and  is  evaluated  at  z=L/2.  Equation  56  can  be  numerically  integrated 
to  give 


[Eq58] 

where  =  ti+]  -  ti,  and  the  load  and  shape  functions  are  a  weighted  average  of  the 
values  at  and  as  follows: 


4)„s(1-a)0(ti)+a4)(t|J 


Pa-(1-a)P(ti)"aP(t,.i)  J 

In  general,  one  should  use  a=0.5,  i.e.,  the  trapezoidal  rule. 


[Eq59] 


It  is  of  interest  to  note  that  for  a  constant  load  P  the  beam  will  creep  linearly  with 
time.  If  it  is  given  that  m=l,  and  that  the  beam  will  creep  to  10  times  its  elastic 
deflection  in  10  years,  then  the  parameter  B  is  given  by  1/12EI  (ksi'-mo.)'^. 


ksi:  kips  per  square  inch;  1  kip  (kilopound)  equals  1,000  lb. 
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Durability  Model  for  Aluminum 

Degradation  of  the  aluminum  beam  will  be  treated  the  same  as  for  the  steel  beam 
except  that  there  will  be  no  surface  protection  model.  Degradation  is  caused  by 
corrosion  of  the  exposed  surfaces,  as  previously  shown  schematically  in  Figure  8.  The 
strength  of  the  aluminum  beam  at  time  t  is  given  by 

Pu(t)  =  I  [K%-2x(t)]  +  4hb[t,-2x(t)]]  [Eq  60] 

where  all  of  the  plate  faces  are  corroded  to  a  depth  of  The  strength  of  the 
aluminum  fy,  the  depth  of  the  beam  h,  and  the  width  of  the  beam  b,  are  taken  as 
constants.  It  is  assumed  that  there  are  no  long-term  deflection  problems  for  an 
aluminum  beam. 

Corrosion  of  Aluminum 

The  aluminum  beam  has  an  I-type  cross-section  with  depth  h  and  flange  width  6.  It 
has  a  flange  thickness  tf  and  a  web  thickness  that,  due  to  the  history  of  temperature 
T,  moisture  0,  and  corroding  chemicals  c,  have  corroded  by  the  amount  x  on  each  face. 
Figure  9,  shown  previously,  illustrates  the  measure  of  material  loss  caused  by 
corrosion.  Corrosion  of  the  aluminum  will  be  modeled  in  exactly  the  same  manner  as 
for  steel  model,  but  the  corrosion  rate  is  much  slower  for  aluminum. 

The  following  diffusion  equation  describes  the  loss  of  plate  thickness  caused  by 
corrosion: 


dx  _ 

^  ’  x2(t) 


[Eq61] 


x(0)  =  0  J 

The  cross-section  starts  in  an  undamaged  condition.  The  corrosion-driving  function 
is  given  by  the  expression 

f,CT.0,c)  =  Djewrn  Mt)r  e  [Eq  62] 


The  parameters  n,  m,  D„,  and  P  govern  the  rate  of  corrosion.  These  parameters  can  be 
established  in  exactly  the  same  manner  as  in  the  model  for  steel  corrosion. 
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The  solution  to  the  rate  equation  for  corrosion  can  be  obtained  by  numerical 
quadrature  to  give  the  loss  of  plate  thickness  at  time  in  terms  of  the  loss  at  time 
t'’i  augmented  by  the  corrosion  that  takes  place  between  those  two  times: 


_i 

x(^.,)  =  [x'(t|)+3Atf,(T„.0„.cJ]® 


[Eq63] 


x(0)  =  0 


Durability  Model  for  FRP 

Creep  Deflection  of  FRP 

FRP  is  subject  to  creep  deflections  caused  by  long-term  loading.  To  model  the 
long-term  creep  behavior  a  Norton  exponential  creep  law  (Kachinov  1986)  was  adopted 
and  assumed  to  apply  to  the  moment-curvature  relationship  of  the  beam  as  follows: 


K(z,t)=iM+BM'"(z,t)  [Eq64] 

where  k  is  the  rate  of  change  of  curvature  of  the  beam,  M  is  the  bending  moment,  and 
El  is  the  flexural  modulus  (the  latter  of  which  may  vary  with  time  due  to  loss  of 
cross-section).  The  parameters  B  and  m  are  the  creep  coefficients. 

For  the  model  problem  M(z)  -  Pz/2,  enabling  the  integration  of  Equation  64  with 
respect  to  the  spatial  dimension  z,  and  resulting  in  an  equation  for  the  rate  of  change 
of  deflection  as  follows: 


W(z,t)  =  (t»(z,t)P(t)+i|)(z)P'"(t) 

where  the  beam  deflected-shape  functions  are  given  by 


<Mz.t)- 


3zL^-4z^ 

48EI(t) 


u)(z)-B(m^2)zL'"-^-2'"-^z'"-" 

(m+1)(m+2)22™*i 


[Eq65] 


[Eq66] 
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The  maximum  deflection  occurs  at  z=LI2.  A  hat  is  used  to  denote  any  variable  that 
both  depends  on  z  and  is  evaluated  at  z=LI2.  Equation  65  can  be  numerically 
integrated  to  give 


w(v,)  =  w(t,)4aIP(Vi)-P(t,)]+iipP:At 


[Eq67] 


where  At  =  -  1,  and  the  load  and  shape  functions  are  a  weighted  average  of  the 

values  at  and  as  follows: 


[Eq68] 


P„=(1-a)Pft)+aPft.d 


In  general,  one  should  use  cc=0.5,  i.e.,  the  trapezoidal  rule. 

It  is  of  interest  to  note  that  for  a  constant  load  P  the  beam  will  creep  linearly  with 
time.  If  it  is  given  that  m=l  and  that  the  beam  will  creep  to  10  times  its  elastic 
deflection  in  10  years,  then  the  parameter  B  is  given  by  1/12E/  (ksi-mo.)'^. 
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Figure  9.  Measure  of  steel  corrosion  depth. 
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4  Modeling  the  Durability  of  a  Cladding 
Element 

Description  of  Modei  Problem 

As  a  model  problem  consider  the  durability  of  a  cladding  element  of  unit  area  exposed 
to  the  environment.  The  surface  material  can  be  paint,  brick,  aluminum,  pol5rvinyl- 
chloride  (PVC),  or  EIFS  (exterior  insulation  and  finish  system).  The  surface  is 
subjected  to  environmental  factors  that  fluctuate  over  time.  Important  among  these 
factors  for  this  model  problem  are; 

•  ambient  temperature  T(t) 

•  ambient  moisture  Q(t)  due  to  humidity  and  precipitation,  expressed  as  a  fraction 
of  saturation 

•  concentration  of  aggressive  chemicals  or  agents  of  decay  c(t),  expressed  as  a 
fraction  of  the  maximum  concentration. 

Serviceability  Criterion 

The  main  concern  for  the  example  cladding  systems  is  serviceability,  which  includes 
both  aesthetic  appeal  and  protection  quality.  The  serviceability  may  degrade  over  time 
because  of  corrosion,  scaling,  cracking,  peeling,  or  other  causes.  Degradation  of 
serviceability  can  be  modeled  with  rate  equations  much  as  was  done  for  the  structural 
element  (see  Chapter  3).  In  each  case,  the  serviceability  will  be  indexed  with  a  single 
parameter  a(t).  The  initial  serviceability  index  will  be  unity,  and  will  degrade  from 
there.  It  is  assumed  that  the  cladding  element  is  serviceable  if 

a(t)  Y  [E*1 

where  y  is  the  serviceability  limit.  One  can  use  this  limit  to  make  repair  decisions  (i.e., 
repair  ifa(t)  <  y). 
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Durability  Model  for  Paint 

The  serviceability  of  the  painted  surface  a(t)  is  governed  by  a  damage  model.  It  is 
assumed  that  the  rate  of  change  of  the  serviceability  index  is  inversely  proportional 
to  the  value  of  the  index,  and  is  driven  by  a  function  of  the  environmental  factors  as 
follows: 


da''_  IpO'.Q.c.v) 
dt  a''!!) 


a''(0  =  gp{‘1v) 


[Eq701 


a^lO)  =  1 


where  the  function  driving  the  change  in  serviceability  is 

fp(T,e,c,v)  =  Dp[T(t)]“[1  +pe(t)l[1  +YC(t)]  [Eq  71] 


and  where 


T(t)^ 


To-T(t) 

T„ 


tEq72] 


is  a  normalized  measure  of  the  temperature.  It  has  been  noted  that  damage  to  paint 
generally  occurs  at  a  greater  rate  for  very  cold  temperatures  and  for  very  hot 
temperatures,  but  slows  for  intermediate  temperatxires.  Here  the  reference 
temperature  is  taken  as  =  60F.  Deterioration  speeds  up  in  the  presence  of  moisture 
and  aggressive  chemicals,  but  proceeds  even  in  a  dry,  chemically  neutral  environment. 
The  constant  is  taken  to  be  0.015  sec  ^ 

The  quality  of  the  surfacing  determines  the  condition  to  which  the  porosity  index 
returns.  The  fimction  used  in  this  model  is  simply  taken  to  be 


gp(q)  =  0.9+0.1q 


[Eq73] 
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The  parameters  a,  p,  and  y  and  are  taken  to  be  functions  of  the  quality  of  the  surfacing 
(i.e.,  material  used,  surface  preparation,  etc.).  The  user  must  specify  the  quality  q  of 
the  initial  surfacing  and  of  any  subsequent  repair  surfacing.  The  quality  q  will  be 
specified  as  a  number  between  0  (poorest  quality)  and  1  (best  quality).  Then  the 
parameters  at  the  vth  repair  can  be  specified  in  accord  with  the  quality  of  the  repair 
using  the  following  relationships: 


Pv  =  1-‘k 

Yv  =  1-‘k 


[Eq74] 


The  maintenance  strategy  then  consists  of  specifying  the  times  and  qualities  of  the 
resurfacing.  For  current  purposes  it  is  assumed  that  the  repair  program  consists  of 
a  single  quality  of  repair  executed  each  time  the  serviceability  index  inequality  is 
violated. 


Equation  70  can  be  numerically  integrated  to  give  the  following  value  of  the  porosity 
at  time  in  terms  of  the  porosity  at  time  after  the  vth  repair: 


a^fCi)  =  /[a''(t,Y  -2fp(T„.0„.c„)At 


a''(C)  =  gp(q) 


[Eq75] 


a“(0)  =  1 


where  At  =  -  t\,  and  the  environmental  factors  are  the  trapezoidal  averages  of 

values  at  and  as  before. 


Durability  Model  for  EiFS 

EIFS  is  similar  to  paint  in  two  important  respects:  it  is  used  as  an  external  cladding 
and  it  is  a  polymeric  material.  Therefore,  for  the  current  work  it  is  assumed  that  the 
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serviceability  of  EIFS  is  essentially  the  same  as  paint,  but  with  different  constants 
used  in  the  degradation  model.  The  constant  is  taken  to  be  0.002  sec'^. 


Durability  Model  for  Aluminum 

Degradation  of  the  aluminum  surface  will  be  caused  by  corrosion,  an  oxidation  process 
the  rate  of  which  increases  exponentially  with  temperature.  The  aluminum  surface 
quality  a(t)  decays  over  time  due  to  the  history  of  temperature  T,  moisture  0,  and 
corroding  chemicals  c.  The  following  equation  describes  the  loss  of  serviceability 
caused  by  corrosion: 


^  =  -f,(T,e,c)a(t) 
at 


[Eq76] 


a(0)  =  1 


The  surface  starts  in  an  undamaged  condition.  The  corrosion-driving  function  is  given 
by  the  expression 


f^(r,e,c)  =  A[0(t)]"[1  +c(t)]"'  e  [Eq  77] 

The  parameters  n,  m,  A,  and  P  govern  the  rate  of  corrosion.  These  parameters  can  be 
established  in  exactly  the  same  manner  as  for  the  steel  corrosion  model.  The  values 
used  in  this  case  were  n=l,  m=3,  A=1.5xl0  ®  sec  \  and  P=450K 

The  solution  to  the  rate  equation  for  corrosion  can  be  obtained  by  numerical 
quadrature  to  give  the  serviceability  at  time 


a(W  =  aft)e 


-fcO’a.ea.Oa)** 


[Eq78] 


a(0)  =  1 
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Durability  Model  for  PVC 

Degradation  of  the  PVC  surface  will  be  caused  by  ultraviolet  radiation  from  sunlight. 
To  simplify  this  discussion  it  is  assumed  that  ultraviolet  exposure  is  a  function  of 
ambient  temperatiu-e.  The  surface  deteriorates  by  an  amount  a(t)  due  to  the  history 
of  temperature  T  and  chemical  exposure  c.  PVC  is  impervious  to  moisture.  The 
following  diffusion  equation  describes  the  loss  of  serviceability  due  to  ultraviolet 
degradation: 


da  ^  fuCiC) 
dt  a(t) 


[Eq79] 


a(0)  =  1 


The  surface  starts  in  an  undamaged  condition.  The  driving  function  is  given  by  the 
expression 


f„(T.0.c)  =  Ae^(')[Uc(t)r 


[EqBO] 


where 


T(t)- 


T(t)-T, 

To 


[Eq81] 


is  a  normalized  measure  of  the  temperature.  Damage  caused  by  ultraviolet  exposure 
generally  occurs  only  at  higher  temperatures  (i.e.,  when  exposure  to  sunlight  is  more 
intensive).  Here  the  reference  temperature  is  taken  as  =  300  °K.  Deterioration  also 
speeds  up  in  the  presence  of  chemicals.  The  parameters  7n=l  and  A  =  0.5x10  ®  sec'^ 
govern  the  rate  of  degradation. 

The  solution  to  the  rate  equation  can  be  obtained  by  numerical  quadrature  to  give  the 
serviceability  at  time 
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[Eq82] 


Durability  Model  for  Brick 


Degradation  of  the  brick  surface  will  take  the  form  of  scaling  caused  by  the  history  of 
temperature  T,  moisture  0,  and  corroding  chemicals  c.  The  following  diffusion 
equation  describes  the  loss  of  serviceability  due  to  scaling: 


$  =  -f.(T.0.c) 

ut 

a(0)  =  1 


[Eq83] 


The  surface  starts  in  an  undamaged  condition.  The  driving  function  is  given  by  the 
expression 


f3(T,e,c)  =  C<T,-T(t)>[e(t)r[1  +c(t)r 


[Eq84] 


where  the  parameters  C,  n,  and  m  affect  the  rate  of  scaling  caused  by  moisture  and 
concentration  of  chemicals,  much  as  in  the  corrosion  model.  The  notation  <x>  means 
that  <:a:>=  x  if  x>0  and  <*:>=  0  if  x<0.  The  value  T^=30S  °K  is  the  fixed  temperature 
above  which  no  freeze-thaw  cycling  is  expected  to  occur.  The  values  used  here  are 
C=1.0xl0'^  sec^K^,  n=3,  and  m=3. 

The  solution  to  the  rate  equation  can  be  obtained  by  numerical  quadrature  to  give  the 
serviceability  at  time  t\^j: 


a(Vi)=a(ti)-2Atfe(T„,e„,c„) 


[Eq85] 


a(0)  =  1 
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5  Durability  Model  Test  Cases 

Test  Case  Overview 

This  chapter  presents  durability  cost  results  of  two  BMDM  test  cases.  Two  types  of 
building  components  were  tested:  a  simply-supported  beam  and  a  cladding  element. 
The  beam  materials  tested  were  reinforced  concrete  and  steel;  the  cladding  materials 
tested  were  paint  and  aluminum.  The  design  criteria  and  durability  models  for  beams 
and  cladding  were  detailed  in  Chapters  3  and  4,  respectively. 

The  material  systems  were  tested  in  two  different  environments  given  the  generic  form 
described  previously  in  terms  of  temperature  T(^),  moisture  content  Qit),  and  chemical 
concentration  c(t): 

T(t)  =  T„+ATsin(nt/6)+6Tti(t) 


0(t)  =  0^+A0sin(7rt/6)+60Ti(t) 


c(t)  =  c„  +  6cTi(t) 


The  two  environments  were  given  the  descriptive  names  “Tropical”  and  “Continental 
U.S.,”  defined  by  the  following  parameters  and  values: 

Tropical 

T„  =  85“F  AT=  10°F  6T=  5'F 

0  =  90%  A0=5%  60=1% 

c„  =  0.2  6c  =  0.1 

Continental  U.S. 

T„  =  40°F  AT=  50'F  6T=  lOT 

0„=  50%  A0=3O%  60=5% 

Cg  =  0.7  6c  =  0.3 
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Figures  12  and  13  show  typical  variations  of  these  environmental  stimuli  over  a  120- 
month  period. 


Results  for  Simply  Supported  Beam 

The  reinforced  concrete  beam  degrades  over  time  through  three  mechanisms: 
reinforcement  bar  (rebar),  concrete  scaling,  and  static  fatigue.  The  beam  is  repaired 
by  restoring  the  member  to  its  original  state.  The  steel  beam  degrades  through  corro¬ 
sion  only;  a  layer  of  paint  is  applied  to  the  beam  to  slow  down  the  corrosion  process; 
but  this  paint  also  degrades  over  time.  The  steel  beam  is  repainted  to  a  user-specified 
quality  whenever  the  paint  (or  surface)  quality  falls  below  a  certain  value.  The  steel 
beam  itself  is  repaired  only  by  replacement. 

Figures  14  and  15  show  some  typical  results  obtained  from  the  BMDM  spreadsheet 
analysis  of  the  reinforced  concrete  beam  subjected  to  both  the  tropical  and  continental 
U.S.  stimuli,  respectively.  It  is  apparent  in  this  model  that  corrosion  of  the  rebar  is 
the  only  mechanism  of  deterioration  present  in  the  tropical  environment,  but  both 
rebar  corrosion  and  concrete  scaling  occur  in  the  continental  U.S  environment.  Table 
1  shows  durability  costs  for  the  beam,  measured  by  taking  averages  of  10  simulations. 
Different  time  periods  (50  and  10  years)  and  discount  rates  (5%  and  20%  annually)  Eire 
considered  for  each  environment.  Over  a  10-year  period,  no  repairs  are  necessary  in 
either  environment,  so  durability  cost  is  the  same  (i.e.,  equal  to  the  initial  cost).  Over 
50  years,  one  repair  is  needed  in  both  environments.  However,  because  the  beam 
subjected  to  the  continental  U.S.  environment  needs  to  be  repaired  sooner  (see  Figures 
14  and  15),  the  durability  cost  for  the  material  is  also  higher.  Note  in  Table  1  that  the 
higher  discount  rate  reduces  the  cost  difference  between  the  two  environments  because 
the  present  value  of  money  at  any  point  in  the  future  decreases  as  the  discount  rate 
increases. 

Figures  16  through  19  show  sample  calculations  for  the  steel  beam.  Figures  16  Emd 
17  show  data  obtained  using  the  tropical  stimuli,  with  two  maintenance  (i.e., 
repainting)  strategies;  low  quality  (0.2,  Figure  16)  and  high  quality  (0.6,  Figure  17). 


Table  1.  Durability  costs  for  reinforced  concrete  beam  (averages  from  10  BMDM  simulations). 


Tropical 

Continental  U.S. 

Cost  ($/m) 

50yrs 

10yrs 

50yrs 

10yrs 

5%  discount  rate 

192.57 

168.44 

205.71 

168.44 

20%  discount  rate 

168.51 

168.44 

168.85 

168.44 
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Figures  18  and  19  show  similar  data  for  the  continental  U.S.  environment.  Note  that 
the  repair  trigger  (i.e.,  the  surface  quality  measure  that  triggers  repainting)  is  set  at 
0.6  in  all  cases.  Table  2  shows  average  durability  costs  for  the  steel  beam  in  the  two 
environments  for  different  maintenance  strategies,  time  periods,  and  discount  rates. 
In  the  tropical  environment,  there  is  little  difference  between  the  costs  associated  with 
a  high-  and  low-quality  maintenance  strategy.  The  low-quality  approach  results  in  one 
repainting  near  the  end  of  the  50  year  time  period,  which  is  inexpensive.  However, 
considerably  more  painting  is  required  in  the  continental  U.S.  if  a  low-quality  mainte¬ 
nance  strategy  is  adopted  instead  of  high-quality  maintenance.  This  results  in 
durability  costs  two  to  four  times  higher  than  those  for  a  high-quality  maintenance 
strategy  in  this  environment. 


Results  for  Cladding 

Paint  and  aluminum  cladding  materials  were  analyzed  by  a  method  similar  to  that 
used  for  the  concrete  and  steel  beams.  As  discussed  in  Chapter  4,  both  t5q)es  of 
cladding  materials  suffer  degradation  of  surface  quality  as  measured  using  a  variable 
between  0  and  1.  Paint  is  maintained  by  repainting  to  a  specified  quality,  and  the 
aluminum  is  repaired  by  replacement. 

Figures  20  through  23  show  the  behavior  of  the  paint  cladding  in  the  two  test  environ¬ 
ments  when  two  repair  strategies  are  used  (low-quality,  0.2,  and  high-quality,  0.8). 
Table  3  shows  average  values  of  the  d\irability  costs  associated  with  these  strategies 
for  50-  and  10-year  periods,  using  discount  rates  of  5%  and  20%.  The  repair  trigger 
was  held  fixed  at  0.6  in  all  cases.  In  the  tropical  environment,  low-quality  repairs  are 
more  expensive  than  high-quality  maintenance  because  no  fimther  repainting  is 
required  in  the  latter  case.  The  opposite  is  true  in  the  continental  U.S.;  in  terms  of 
durability  cost  the  relatively  high  frequency  of  low-quality  repairs  is  offset  by  the  low 
per-repair  cost  as  compared  to  high-quality  (and  high-cost)  repairs  which  are  also 
needed  frequently.  Thus,  according  to  BMDM  as  run  for  the  continental  U.S.  environ¬ 
ment,  frequent  cheap  repairs  would  be  preferred  to  fewer,  more  expensive  repaintings. 

The  behavior  and  durability  costs  associated  with  the  aluminum  cladding  are  shown 
in  Figures  24  and  25,  and  in  Table  4.  The  repair  trigger  in  this  case  also  is  fixed  at  0.6. 
Two  repairs  are  required  in  a  50-year  period  in  the  continental  U.S.  environment, 
compared  to  one  repair  required  in  the  tropics.  Therefore,  the  durability  cost  is  higher 
in  the  continental  U.S.  over  the  50-year  period.  The  durability  costs  are  equal  to  the 
initial  cost  over  a  10-year  period  in  all  cases. 
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Table  2.  Durability  costs  for  steel  beam  (averages  from  10  BMDM  simulations). 


Tropical 

Continental  U.S, 

SOyrs 

lOyrs 

SOyrs 

10yrs 

Cost  ($/m) 

0.2  main 

0.6  main 

0.2  main 

0.6  main 

0.2  main 

0.6  main 

0.2  main 

0.6  main 

5% 

discount 

rate 

50.35 

49.53 

49.53 

49.53 

206.23 

50.26 

115,79 

49.53 

20% 

discount 

rate 

49.54 

49.53 

49.53 

49.53 

90.24 

49.54 

84.41 

49.53 

Table  3.  Durability  costs  for  paint  cladding  (averages  from  10  BMDM  simulations). 


Tropical 

Continental  U.S. 

50  yrs 

10  yrs 

50  yrs 

10  yrs 

Cost  ($/m) 

0.2  main 

0.8  main 

0.2  main 

0.8  main 

0.2  main 

0.8  main 

0.2  main 

0.8  main 

5% 

discount 

rate 

18.57 

10.00 

10.00 

29.97 

64.62 

17.66 

32.78 

2o% 

discount 

rate 

11.52 

10.00 

11.17 

10.00 

14.50 

21.98 

13.68 

19.99 

Table  4.  Durability  costs  for  aluminum  cladding  (averages  from  10  BMDM  simulations). 


Tropical 

Continental  U.S. 

Cost  ($/m) 

50  yrs 

10  yrs 

50  yrs 

10  yrs 

5%  discount  rate 

71.28 

50.00 

97.4 

50.00 

20%  discount  rate 

50.22 

50.00 

51.57 

50.00 
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Discussion 

In  general,  different  material  behavior  and  durability  costs  were  observed  for  the  two 
test  environments.  The  concrete  beam  simulations  demonstrated  that  the  nature  of 
the  damage  can  depend  on  the  environmental  stimuli.  Concrete  scaling  was  observed 
only  when  the  beam  was  subjected  to  the  continental  U.S.  stimuli,  but  rebar  corrosion 
was  present  in  both  cases. 

It  was  also  observed  that  the  selection  of  the  most  economical  maintenance  strategy 
can  depend  on  the  environment.  The  BMDM  simulations  suggest  high-quality  mainte¬ 
nance  of  the  steel  beam  for  the  continental  U.S.  stimuli,  but  on  the  other  hand,  the 
choice  of  maintenance  strategy  in  the  tropical  environment  does  not  produce  any  no¬ 
ticeable  cost  differences.  According  to  BMDM  simulations  for  paint  cladding,  high- 
quality  maintenance  is  appropriate  in  the  tropics,  but  low-quality  repainting  is  the 
most  economical  maintenance  strategy  for  the  continental  U.S.  environment. 

The  simulations  and  observations  at  this  point  in  the  research  confirm  that  BMDM 
can  qualitatively  predict  and  compare  the  durability  cost  of  various  building  materials. 
The  algorithms  upon  which  BMDM  is  based — ^formulas  representing  materials 
durability  and  environmental  stimuli — ^would  need  to  be  extended  and  refined  before 
BMDM  can  quantitatively  model  complex  materials  degradation  and  environmetal 
processes.  However,  the  findings  of  this  study  clearly  indicate  that  BMDM  could  be 
developed  into  a  finished  tool  to  help  facility  designers  and  engineers  make  better  deci¬ 
sions  in  selecting  building  materials  for  optimal  life-cycle  economy  and  serviceability. 
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Figure  13.  Typical  continental  U.S.  environmental  stimuli. 


Figure  14.  Behavior  of  concrete  beam  subjected  to  tropical  environmental  stimuli. 
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Figure  15.  Behavior  of  concrete  beam  subjected  to  continentai  U.S.  environmentai  stimuii. 


Figure  16.  Behavior  of  steei  beam  subjected  to  tropical  stimuli  (low-maintenance  strategy). 


steel  Beam 
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Figure  17.  Behavior  of  steel  beam  subjected  to  tropical  stimuli  (high-maintenance  strategy). 


steel  Beam 
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Figure  18.  Behavior  of  steel  beam  subjected  to  continental  U.S.  stimuli  (low-maintenance  strategy). 
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Figure  19.  Behavior  of  the  steel  beam  subjected  to  continental  U.S.  stimuli  (high-maintenance  stragety). 


Figure  20.  Behavior  of  paint  cladding  subjected  to  tropical  stimuli  (low-maintenance  strategy). 


Figure  21 .  Behavior  of  paint  cladding  subjected  to  tropical  stimuli  (high-maintenance  strategy). 


Figure  22.  Behavior  of  paint  claedKig  subjected  to  continental  U.S.  stimuli  (low-maintenance  strategy). 
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Figure  23.  Behavior  of  paint  cladding  subjected  to  continental  U.S.  stimuli  (high-maintenance  strategy). 
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Aluminum  Cladding 


Figure  25.  Behavior  of  aluminum  cladding  subjected  to  continental  U.S.  stimuli. 
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6  Conclusions  and  Recommendations 


Summary 

A  Building  Materials  Durability  Model  (BMDM)  has  been  developed  and  implemented 
as  a  computer  simulation  to  calculate  the  net  present  value  of  a  building  component 
as  estimated  from  the  component’s  first  cost  and  the  estimated  cost  of  required  future 
repairs.  BMDM  accepts  arbitrary  inputs  of  environmental  stimuli,  determines  the 
durability  for  a  token  system  having  the  nominal  design  requirements  of  the  real 
structural  component  in  service,  considers  simultaneous  accumulation  of  damage  due 
to  multiple  environmental  phenomena,  predicts  the  time  to  dysfunction  or  failme,  and 
estimates  total  durability  cost  as  the  sum  of  first  costs  and  the  current  cost  of  future 
repairs. 

The  simulations  discussed  in  Chapter  5  show  that  BMDM  can  be  used  to: 

•  compare  the  life-cycle  costs  of  two  candidate  materials  for  a  given  application 

•  assess  the  consequences  of  different  repair  strategies 

•  compare  durability  of  a  material  in  different  environments  (locations) 

•  choose  the  best  diuability  strategy  for  short-  and  long-term  applications 

•  consider  the  time  cost  of  money  in  making  durability  strategy  decisions. 

It  is  concluded  that  BMDM  could  be  used  as  the  basis  of  a  Monte  Carlo  simulation  of 
a  material’s  total  durability  cost  and  its  uncertainty. 

Work  to  date  represents  an  initial  development  attempt  and  proof  of  concept. 


Recommendations 

To  fulfill  its  potential  as  a  materials  selection  and  costing  tool  for  facility  designers, 
BMDM  requires  refinements  and  extensions.  It  is  recommended  that  future  work  on 
BMDM  include  the  following: 

•  improving  models  for  a  given  material  and  application  ixsing  careful  comparisons 
with  observed  behavior 
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•  considering  additional  token  systems 

•  surveying  and  cataloguing  the  durability  problems  for  construction  materials 

•  building  up  a  library  of  damage  models  for  needed  applications 

•  increasing  the  sophistication  of  the  environmental  inputs  and  consider  in  more 
detail  the  relationship  between  measured  environmental  factors  and  the  actual 
variation  in  these  quantities  experienced  by  the  material 

•  altering  the  current  model  to  permit  the  initial  quality  of  the  building  material 
or  component  to  be  a  random  variable 

•  generalizing  the  framework  for  flexibility  with  respect  to  geometry  and  loading 

•  considering  transport  models  (e.g.,  for  moisture,  chemicals  and  elements,  heat, 
etc.)  in  more  detail. 
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